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Abstract. The origin of large-scale magnetic fields in cosmic structures and the intergalactic 
medium is still poorly understood. We explore the effects of non-minimal couplings of elec- 
tromagnetism on the cosmological evolution of currents and magnetic fields. In this context, 
we revisit the mildly non-linear plasma dynamics around recombination that are known to 
generate weak magnetic fields. We use the covariant approach to obtain a fully general and 
non-linear evolution equation for the plasma currents and derive a generalised Ohm's law 
valid on large scales as well as in the presence of non-minimal couplings to cosmological 
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vational bounds on such couplings, we conclude that modifications of the standard (adiabatic) 
evolution of magnetic fields are severely limited in these scenarios. 
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1 Introduction 

Our Universe is apparently magnetised on virtually all different length-scales probed by 
astronomical observations. The strength of magnetic fields in galaxies and galaxy clusters 
is of the order of /iGauss [1—4], and a lower bound of the order of 10 -15 Gauss on coherent 
magnetic fields in the intergalactic medium has also been reported [5-9] . The origin of these 
large-scale magnetic fields still remains unclear. Speculative mechanisms, active at all stages 
of the evolution of the Universe, have been proposed which often require the breaking of the 
conformal invariance of electromagnetism (EM). For recent reviews see refs. [10-13]. 

Besides these exotic possibilities, it is known that weak magnetic fields are generated 
around photon decoupling. In the pre-recombination era, the vorticities of electrons and 
protons evolve slightly differently yielding a net circular current that generates magnetic 
fields, an effect known as the Harrison mechanism [14, 15]. However, no magnetic fields are 
produced at first order in cosmological perturbations, even in the presence of active sources 
of vector perturbations [16]. At second order in both perturbations and the tight coupling 
approximation of the photon-baryon interactions, the mechanism does generate magnetic 
fields of about 10~ 29 Gauss on Mpc scales [17-26]. 

In this paper, we investigate whether such seeds could be boosted at late times due to 
non-minimal (effective) interactions of EM. In particular, we allow for different extensions 
of EM to embrace couplings with a scalar or pseudo-scalar field, for instance, dark energy 
or a cosmological axion. In this context, we review the hydrodynamical evolution of the 
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primordial plasma from electron-positron annihilation until today by carefully taking into 
account the effects of a possible non-minimal theory. 

First of all, in section 2 we provide a brief review of the plasma processes occurring 
from the pre-recombination era to today by looking at the evolution of currents, relevant 
on cosmological scales at second order. Then, we derive a generalised Ohm's law that is 
also valid in the context of non-minimally coupled EM. By inspecting the different times- 
cales involved in the current evolution we briefly discuss its validity on very large (sub- and 
super- horizon) cosmological scales. This section should serve as a reference for cosmological 
plasma conductivity formulae and physics. In section 3, we introduce the generic form of 
the Lagrangian allowing for different non-minimal couplings and discuss their effects on the 
magnetic field evolution at large scales. Finally, in section 4, we summarise our results and 
discuss the conclusions. Appendix A lists our conventions, and appendix B contains the 
details of the 1+3 covariant approach for the description of the fluid and current evolution 
in the plasma. 

2 Generalised Ohm's law for non-minimal electromagnetism 

Circular currents induce magnetic fields. In standard EM this can easily be seen by taking 
the curl of Ampere's circuital law and substituting the electric field using Faraday's law of 
induction. One arrives at a wave equation for the magnetic fields, OB = — V x J, that 
is sourced by the curl of the current. In curved space and, in addition, in scenarios with 
non-minimal couplings of EM the situation is more complicated. However, the curl of the 
standard electric current i.e. the net flux of charged standard model fermions will still act 
as a source in much the same way. 

To solve for the evolution of the electromagnetic fields and the plasma one typically 
applies Ohm's law. It represents a shortcut to solving Maxwell equations by providing a 
link between the current and the electromagnetic fields, J = a(E + vxB), where v is the 
bulk velocity of the plasma and a is the conductivity that has units of an inverse length. 
Ohm's law derives from the evolution equation of the current and reduces to this simple form 
after assessing the relevant evolution and interaction timescales. For a general relativistic 
derivation of the current evolution and a generalised Ohm's law for certain situations (see 
ref. [27] and references therein). 

In this section we carefully review the evolution of currents and derive a generalised 
Ohm's law relevant on cosmological scales and also valid in the context of non-minimally 
coupled EM. We follow the derivation of the current evolution equation given in ref. [27], and 
assess the different interaction terms similarly to refs. [16, 21, 22, 24]. 

2.1 Relevant particle species and their interactions 

We aim at deriving the evolution equation of electrical currents on cosmological scales valid 
from temperatures T ~ 500 keV (or z ~ 10 9 ) corresponding to electron-positron annihilation 
all the way down to today. Next to the classical electromagnetic fields described by the 
Faraday tensor (F), the constituents of the plasma that we need to follow are the photons 
(r), the free electrons (e), the free protons (p) and the neutral atoms such as hydrogen and 
helium (H), that start forming due to recombination. In this discussion we can safely ignore 
the neutrinos and dark matter. However, ultimately we are only interested in the dynamics of 
electric currents from electrons and protons, and the evolution of the totality of the massive 
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standard model particles, referred to as baryonic matter (label b), not to be confused with 
actual baryons in the particle physics sense (£>). 

The photons interact through Thomson scattering with the free electrons and protons, 
Crs, where s = e, p. The electrons and protons interact in addition through Coulomb scatter- 
ing, Ce P , a.k.a. Rutherford scattering. Finally, the free charges feel an electromagnetic force, 
Cg F , while the neutral atoms only interact gravitationally: 

V a Tf* = -C? p + C$ e +C; F , (2.2) 

v a rr = -c&-c£e + C> ( 2 - 3 ) 

V Q T£ a = 0. (2.4) 

The conservation equation for the baryonic matter reads 

V Q 7T = v « ( T r + T£ a + Tg*) = -C& - C& + C; F + e eF , (2.5) 

and the Coulomb scattering drops out. 

The energy-momentum tensor of the electromagnetic fields derives from its potentially 
non-minimal action, see next section. In case of the standard Maxwell theory, X7 a T F a = 
—F^ a j a , with the 4-current j^. In this section, we consider a Lagrangian with a canonically 
normalised kinetic term, — (1/4)F 2 , but take into account possible couplings to some scalar 
fields. Effectively, this means we can write 

V a T£ a = W - F» a j Q , (2.6) 

where ZC M combines all contributions to the conservation equation due to the interactions with 
fields beyond the standard model. In addition, we need to allow for a rescaled coupling to 
the standard model fields, a — > a = a(l + 5a), where 8a depends on the coupled fields. The 
stringent bounds on the fine-structure constant, a, translate into the requirement 5a -C 1 and 
very little space-time dependence, such that we may neglect derivatives of 5a. The modified 
a also enters the Thomson and Coulomb collision cross sections. Clearly, the interaction 

terms must add up to zero, such that ^ s V a T^ a = 0, if we also include all fields coupled to 

piiv 

Let us briefly recapitulate how the energy-momentum transfer terms come about in 
the above conservation equations. In kinetic theory, the photons are described by their 
phase-space distribution function that follows a Boltzmann equation with collision terms due 
to the scattering of photons into and out of a given phase-space element. The Boltzmann 
equation is expanded into a multipole hierarchy, where the multipoles are integrations of the 
distribution function over momentum, see e.g. [28]. The first two multipole equations are the 
energy-momentum conservation equations (2.1). They take the usual hydrodynamical form. 
The energy density p r , the energy flux qfi , and the anisotropic stress tt^ of the photons can be 
interpreted as effective fluid variables that describe the macroscopic, collective properties of 
the photons, as opposed to the distribution function that describes the microscopic properties 
of individual photons. The interaction terms, Cre and Cr P , reflect the energy-momentum 
transfer due to the collection of local collisions of the photons with the electrons and protons, 
respectively. In ref. [28] they are computed up to second order, see also ref. [29]. They are 
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split relative to the frame as Cr S = C rs u^ + and for the electrons it is found that 



C re = n e (JT 



n e GT 



p r v\ -v ea q" + 0(3) 



(2.7) 
(2.8) 



where ot = 87ra 2 /(3m 2 ) is the Thomson cross-section. The relative velocity Ve stems from 
the decomposition of the rest frame 4-velocity w.r.t. the generic frame: ut = J e ( u ^ + v e ) 
with UpVe = and j e = (l — u 2 ) -1 / 2 , see appendix B for details. For the protons we have 



nil = al nil 



P 



(2.9) 



Here the factor /3 2 ~ 3 x 10 7 makes up for the mass-dependence of the cross-section and 
strongly suppresses the proton Thomson term as compared to the electron one. Notice that 
the momentum transfer rate, Cre , is first plus second order in cosmological perturbations 
and also includes perturbations of the number density n e . The energy transfer rate, C re , 
starts only at second order in the relative velocity. 

The non-relativistic massive particles are well described by a hydrodynamical treatment 
with vanishing heat flux, viscosity and anisotropic stress. In other words, their Boltzmann 
hierarchy truncates after the dipole, and the Coulomb collision term is derived and interpreted 
in the same way as the Thomson terms discussed above. Following refs. [29] and [30] we find 
for non-relativistic electron-proton collisions 



On 



C$> = n e n p acVTem e (1 + fi) 



< + 0(3) 



(2.10) 
(2.11) 



with the Coulomb cross section oc = 47rA ep a 2 /T 2 and the thermal electron velocity VTe = 
{T e /m e ) l l 2 . The Coulomb logarithm is X ep = 23 + \ ln(T e 3 /n e ) for the temperatures of 
interest [30] . We find A ep ~ 34 at early times before it slowly falls off when the baryons start 
cooling adiabatically, see below. Again, the momentum transfer is first plus second order in 
perturbations including density perturbations, while the energy transfer only starts at second 
order in the velocities. 

Finally, the electromagnetic force terms, C^ F , are not actually due to collisions, but can 
be consistently included as force terms in the Boltzmann equation [29]. The 4-current of 
standard model charges, here the free electrons and protons, reads 



r 



([■ 



n,, 



n e \vr + n p v f 



n e v£) , 



(2.12) 



see appendix B.4 and particularly eq. (B.18). Even though allowing for non-minimal EM, 
we take the interaction between standard model charges and the electromagnetic fields to be 
canonically normalised, as discussed below eq. (2.6). This means we can write the electro- 
magnetic force terms in their standard form 



C s f 



e s n s vfE a , 



C 



</■<} 

sF 



Ef + e^VsaBt 



(2.13) 
(2.14) 
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not forgetting that e s = ±e, the charge of particles of type s, may depend slightly on 
the fields that interact non-minimally with F^ v . But again, we can neglect derivatives of 
e = \J 47ra(l + 5a). The electric and magnetic fields measured in the generic frame are 
projected from the Faraday tensor as follows 

E ll = F IM u a , B„ = ^ a pF a P, F flu = 2u^E^+e^ a B a . (2.15) 

The evolution of the electric and magnetic fields is governed by the (potentially modified) 
Maxwell equations, but we will not use them for the moment. 

Before moving on to the current evolution equation, let us summarise the homogeneous 
and isotropic thermal evolution of the plasma. The formation of neutral hydrogen and helium 
due to the protons capturing electrons is followed in terms of the ionisation fraction x e . It 
relates the average number density of the free electrons to the average baryon number density, 

n e = x e ns , fiB = n P + uh = n e + fin . (2-16) 

The Universe is charge neutral on average, so we have set n e = n p . The evolution of x e starts 
off at 1 at early times, and drops rapidly when the recombination processes are falling out 
of equilibrium, roughly at T rec ~ 0.3 eV. When recombination is nearly completed, around 
z ~ 1010, it settles at 



1.2 x KT 5 , (2.17) 

il b n 

until reionisation of the plasma due to the radiation from the first stars increases x e back up 
to about 0.1, starting roughly from z ~ 10. We employ a numerical fit to the evolution of x e 
given in eq. (3.201) of ref. [31]. The baryon number density, ng, derives from the constant 
baryon to photon ratio, ns/nr = 2.6 x 10 — 8 0{,/i 2 , and the well known result n r = 2£(3)a -3 . 
The average photon energy density is p r = (ir 2 /15)a -4 . Finally, the non-relativistic species 
have negligible average pressure and an energy density of p s — m s n s . 

To compute the Coulomb collision rate we also need to know the evolution of the baryon 
temperature (incl. electrons). As long as Thomson scattering is efficient the baryons are in 
good thermal equilibrium with the radiation. Then we have T b ~ T r = T cm b/a with the 
CMB temperature today T cm b ~ 2.726 K. However, after sufficient dilution of the plasma, 
the electrons do not meet enough photons anymore to keep thermal contact, and the baryons 
start to cool adiabatically. When this happens can be seen from the evolution equation for 
the baryon temperature 

f b + 2HT b = l^^a T n e {T r -T b ) , (2.18) 
3 m e p b 

with the mean particle mass p. This evolution equation is derived from the first law of 
thermodynamics with Thomson scattering as a means of transferring heat from the photons 
to the baryons, see for instance eq. (69) of ref. [32]. With p b ~ pn b = p(l + x e )nB and 
n e = XefiB we can see that it is the rate 

r b = -^- 2.19 

3 m e 1 + x e 

that governs the thermal evolution of the baryonic matter. It drops below the Hubble rate 
at a redshift of about z ~ 120, when the baryons decouple from the photons and T b oc I /a 2 . 
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2.2 Evolution of currents 



Let us now write and assess the evolution equation for the spatial part of the total current, 



J'' 



= haj a , with the 4-current given in eq. (2.12). In appendices B.3 and B.4 we 
give details on the derivation. From the energy and momentum conservation equations we 
derive an evolution equation for the bulk velocity of each species (B.16). These are then 
used together with the evolution of the number densities (B.20) to arrive at the fully general 
non- linear evolution equation for the total spatial current (B.21). This equation and its 
derivation was also given in ref. [27], but with a minor confusion about the pressure of the 
plasma components. 

We spell out the equation for the electron-proton system, and insert the electromagnetic 
force and the elastic Thomson and Coulomb collision terms from above. Even though the 
electrons and protons are non-relativistic, we do not a priori neglect their pressure, but 
rather realise that it only arises at second order in perturbations, p s ~ n s T s <C p s , for 
s = p , e. We replace the pressure perturbations by means of the electron and proton adiabatic 
sound speeds, c e and c p , see appendix B.5. For the energy density, we use p s = m s n s + 
0(2). Finally, keeping the electromagnetic fields non-perturbative, but neglecting velocity 
and density perturbations of third and higher order, we arrive at 



-eacVTei^ + PY 



J a -V a {vp£ + v«J£)-e{ 



n n — n e ) W 



^pJe ^eJn 



4(Jt 

3m P L' 



n\ n 2 p 



Pp 



v e J e j3vip Jp 



E, 



a £ 

m, 



fia/3 



Jr. 



p a 



Br- 



(e + v Q ^)j£ + v^ pa j p Q ) 



-eripC p + en e c e h C(o) . 



{e + v a v«)j? + w{v ea j?) 



Pe 



(2.20) 



Here the partial electron and proton currents are Je 
The effective Thomson photon current is defined as 



-en e Ve and Jp 



3 

Ap r 



e(n e - Pn p )q? - (J? + (3 3 J p Q ) i$ , 



respectively. 



(2.211 



where qfi and itr v are the photon energy flux and anisotropic stress, respectively. 

The general relativistic current evolution equation (2.20) includes all effects that arise 
up to second order in perturbations, also taking number density fluctuations into account. 
The first line shows all geometric and kinetic effects, comprising the isotropic, shear and 
vortical expansion, non-linear plasma inhomogeneities, and the acceleration of local charge 
fluctuations. The second line describes the effect of Coulomb and Thomson scattering on 
the total current. The third line shows how electric fields generate currents, and it accounts 
for the Hall effect, a coupling between the current and the magnetic field. Finally, the 
last two lines in eq. (2.20) are due to the pressure perturbations (and their transformation 
into the generic frame that we employ to describe the current and the electromagnetic 
fields) . In most treatments of magnetic field generation in the recombination era the pressure 
perturbations have been neglected [18, 21, 24], even though they are known to generate 
magnetic fields on small scales, an effect known as Biermann battery, see e.g. [12]. 
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To better understand the evolution of the current, let us rewrite and simplify eq. (2.20) 



a bit. For this purpose we define the average number density n 



n e 



n. 



x e riB, and 



introduce the Coulomb rate, Tc, the Thomson rate, IV, and the plasma frequency, w. 



c = ena c v Te 



3m P 



uj p = (e 2 n/m e ) 



1/2 



(2.22) 



We change variables for the partial currents and number densities by defining the centre of 
mass velocity and number density fluctuation 



u M = 

u pe — 



4 + Pv'e 

1 + /3 



n(l + /3) 



1 . 



(2.23) 



It is often assumed that the centre of mass velocity and density coincide with the velocity and 
density of the baryonic matter. However, we point out that this is not necessarily correct at 
second order and needs to be checked (numerically) by treating the electrons, protons, and 
neutral atoms separately to compare Vp e with Vjj and v£, etc. Next we define the centre of 
mass current and the local charge number fluctuation 



TV 



en (1 + S pe ) v^ e , 



A = 



n. 



n 



(2.24) 



Using these definitions, we replace the velocities, the partial currents and number densities 
in eq. (2.20) and keep terms up to second order in perturbations. The terms that come from 
the pressure perturbations are strongly suppressed even though c 2 QJg and c^V M /9 s come in 
linearly. They are effectively second order contributions because the sound speed is so small, 
Cg = p e /p ~ (5/3) Tb/m e and c 2 = f3c 2 from eq. (2.18). Therefore, we can safely neglect the 
terms c 2 [ (V a vf) Jg + V M (u SQ ,J°)] as compared to the similar terms in the first line of eq. 
(2.20). Next, to facilitate the qualitative discussion we also consider the electric and magnetic 
fields as being small fluctuations. We neglect terms 0(v 2 E : A 2 E, AvB), even though these 
may have to be kept for a fully consistent treatment [33]. Finally, we neglect (3 = m e /m p as 
compared to 1. With these substitutions and simplifications we arrive at 



•i (4 + 3^)0 J" 



2 



[l + 6 pe )J»-AJ£ e 



en \ pe J 
(1 + S r - p5 pe )J» - (1 + 6 r )J£ + J<? 



1 + 5, 



pe 



A 



P £ pa/3 



en 



pe a 



+c 2 e ejp 1 e + enc 2 e V» [5 pe -A)+ C(3) . 



(2.25) 



where 5 r = p r /p r — 1 is the photon density contrast. In this approximation the effective 
photon current reads 



4p r 



en(l + S pe - A)< + ( J p Q e - J a ) n? a 



(2.26) 



In eq. (2.25) we can see that both Coulomb and Thomson scattering damp currents as well as 
drive them through the terms oc Jp e and Jr. Pressure gradients and electric fields generate 
currents and magnetic fields curl them. 
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We remind the reader that the only impact of non-minimal couplings of EM fields 
on the evolution of standard model currents is an effective space-time dependence of the 
fine-structure constant and therefore the electric charge. We did not use the Maxwell field 
equations, which is the reason why non-standard currents do not appear in the current 
evolution equation. In fact, we did not use the Einstein field equations either, and therefore 
these results also hold in scenarios with modifications to gravity, such as the presence of 
couplings of EM to curvature invariants. 

2.3 Generalised Ohm's law 

Let us now estimate the different timescales involved in the evolution of currents. The main 
features of eq. (2.25) remain if we linearise it in a FLRW background. The proton-electron 
centre of mass current, Jp e , is only present because we worked in an arbitrary frame, Let 
us now choose to be the centre of mass frame, which is very close to the baryon frame, 
then Jp e = 0. We express the spatially projected 4- vectors with respect to a comoving basis, 
i.e. = a -1 (0, J) etc. Then we find the following linear current evolution equation 

j + (4H + Tc + TT)J = -T T J r + <4 E > ( 2 - 27 ) 

where the Hubble parameter H comes from O = 3H + 0(1). To write a generalised cosmo- 
logical Ohm's law in a familiar form, we approximate the time derivative of the current by a 
characteristic timescale of the problem, r, through J ~ t _1 J. For large-scale fluctuations of 
physical correlation length L larger than the Silk damping scale, typically r ~ min(L, H -1 ). 
Then we can write the linear cosmological Ohm's law as 

(Vr + 4 VH + vc + Vt)J + VtJt ^ E , (2.28) 
where the various resistivities are 

Vr = —r > VH = ~2 ) VC = —2, VT = 2 i ( 2 - 29 ) 

w p w p w p 

and have the dimensions of time. The resistivities quantify the efficiency of generating cur- 
rents by the presence of an electric field. We compare the resistivities as functions of the scale 
factor in the upper panel of figure 1 and notice that rjx 3> r]c for a < 3 x 10~ 6 , well before 
recombination, and then r\T <C r]c until today. The resistivities due to the expansion, r/n, 
and the characteristic evolution timescale, rj T , are always negligible, i.e. the time derivative 
in the current evolution equation can safely be neglected. We can write eq. (2.28) in the form 

J + (JEVT Jr ^ <T E E , (2.30) 

where oe = (vc + ^t) -1 is the electric conductivity of the plasma with the dimensions of an 
inverse length. 

On the other hand, we should regard the effective photon current, J r , as the external 
source, rather than the electric field. We need to compare J and &eE to see which one is 
generated by the photon current. To do so we have to use the Maxwell equations. Let us 
assume standard EM for the rest of this discussion. Because of conformal invariance we can 
rescale the electromagnetic fields and the currents according to 

E = a 2 E , B = a 2 B , J = a 3 J, J r = a 3 J r , (2.31) 
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10~ 8 10" 6 10" 4 0.01 1 



a 

Figure 1. Upper panel: we compare the various resistivities, defined in eq. (2.29) for a comoving 
correlation length of 1 Mpc, that are relevant in the cosmological Ohm's law (2.28). Notice that 
i]t S> f]c f° r i^3x 10~ 6 , well before recombination, and then t/t <C 77c until today. Lower panel: 
we plot the comoving conductivity a = aj(r\c + t)t)- 

and use the flat space-time Maxwell equations for the twiddled vectors: V-B = 0, V-E = a 3 Q, 
Vx B = E' + J, and Vx E = —B'. The prime denotes conformal time derivative. From the 
curl of Ampere's law together with the induction equation we find the sourced wave equation 
B" — V 2 B = Vx J. Finally, we take the curl of eq. (2.30) and use the sourced wave equation 
to find 

B" - V 2 B + ctB' ~ -— V x J r , (2.32) 

a 

where the comoving conductivity is a = cute, see figure 1. On comoving scales larger than 
(T _1 one concludes that the first two terms on the left hand side can be neglected w.r.t. 
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the third term. This is equivalent to dropping J with respect to c?eE in eq. (2.30). In 
this case eq. (2.32) can be integrated to compute the magnetic field, which was done in ref. 
[24] . They find the photon current to be an active source through recombination until about 
z ~ 300, even though Thomson scattering starts to become inefficient already well before 
recombination, see figure 1. Magnetic fields of ~ 10 -24 Gauss on large scales are found, see 
the beginning of next section. 

Finally, one may ask whether the cosmological Ohm's law is applicable on very large 
and super- horizon scales. The resistivities that relate the current to the electric field are due 
to scattering processes. Thus, one might conclude that the resistivities are not well defined 
on super-horizon scales due to causality - particles in different Hubble volumes are never in 
causal contact with each other. However, let us emphasise that the general current evolu- 
tion equation (2.20) was derived from the first two multipole equations of the Boltzmann 
hierarchies for the photons, electrons and protons. The collision and EM force terms res- 
ult from passing from the microscopic kinetic theory to the macroscopic fluid description 
of the plasma, through integration over momentum space. They are local, space-time de- 
pendent terms that describe the local energy-momentum transfer rates for the effective fluid 
description, and not interactions between spatially separated particles. Let us assume that 
inflation imprinted super-horizon correlations into the photon, electron and proton distri- 
butions. Then these correlations are passed through momentum transfer to the EM fields, 
even on super-horizon scales, as causality is built into the relativistic Boltzmann hierarchy. 
In other words, inflation ensures that the charge flux in two disconnected Hubble regions is 
correlated intrinsically. Therefore, the EM fields generated by the correlated currents will be 
correlated on super-Hubble scales as well. Consequently, the resistivities, the conductivity 
and Ohm's law are still well defined on super-horizon scales. 

Notice that in the case of non-minimal EM the form of the generalised Ohm's law 
would not change, as we have not used the Maxwell field equations to derive it. Only, 
the fine-structure constant that goes into the resistivities would become slightly space-time 
dependent, as we considered a canonically normalised theory. 

3 Scenarios of modified electromagnetism 

The original motivation for looking at the possibility of late-time couplings between EM and 
cosmological (pseudo-)scalars came from explaining the observed cosmological magnetisation. 
Some time after the completion of recombination the Universe is populated by seed fields 
whose stochastic distribution is characterised by a power spectrum of the form [24] 



where k eq = % eq ~ 0.01/i/Mpc is the mode that enters the horizon at matter-radiation 
equality. This spectral shape is reached roughly at redshift z ~ 300, when the non-linear 
processes saturate and the magnetic fields continue their evolution through adiabatic ex- 
pansion, B tx a -2 , in the absence of any other couplings and effects. At a = Wa eq , where 
a eq = Q r /£l m « 1/3240, the physical field strength on a scale k = 20k eq is found to be 
~ 10 -24 Gauss, see right panel of figure 4 in ref. [24]. Clearly, these fields are not going to 
be strong enough if they are to explain observations. One can then think of adding cosmo- 
logical (pseudo-)scalar fields to this picture, whose interactions with EM might act only at 
late times and perhaps boost these weak seeds to more reasonable amplitudes. What we will 




for k <C k, 
for k » k, 



eq 



eq 



(3.1) 
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find is that, even in the emptiest Universe, the very low resistivity of the plasma will prevent 
meaningful modifications of Maxwell theory from producing any effect. 

We thus consider a scalar or pseudo-scalar field <\> (be it the quintessence, dilaton, axion 
or any other effective cosmological field) with generic Lagrangian density = V(X, </>) and 
kinetic term X = — ^d a 4>d a (j) that is minimally coupled to gravity but non-minimally to EM, 
as given in the following general action [10, 34, 35] 



S[<P,A^]= / d 4 ^ 



^R + V{x,4>)-\f\4>)F a ^ 

+ ~g(<t>)F a pF a P - ^m 2 ^)A a A a +j a A c 



(3.2) 



Here F^ = d fl A v — d u A^ is the EM field tensor and F^ u is its dual, defined as F^ u = 
\rf ,va P F a p with r]^ ua P being the totally antisymmetric Levi-Civita tensor with rj 0123 = 
+ |g| -1 / 2 . The coupling functions f(4>) and g((j>) are dimensionless while m(<f)) represents 
an effective mass for the photon. 

Since the canonical normalisation of the kinetic term of the vector potential leads to a 
redefinition of the electron charge as e — > e/f, the term f 2 {4>) is directly responsible for the 
spatiotemporal variation of the fine-structure constant a. The axial coupling function g((j)) 
arises naturally in theories where a pseudo-scalar is present, such as the axion. It is modeled 
on the analogue anomalous coupling between photons and the axial neutral mesons of QCD. 
Note that we only consider the field coupling to EM and not the derivative couplings as they 
will presumably be suppressed by the scale associated with the derivative interaction. 

The Euler-Lagrange equations in compact form are 

fd a F ail + 2fd a fF apL + 8 a gF afl - m 2 A^ + j„ = . (3.3) 

For the current, we assume a simple Ohm's law, = oeF^ = <te(0, — A^/a), and define 
the comoving conductivity again by a = clge- Here we neglect the effective photon current 
discussed in section 2.3, as it is irrelevant only for z < 300. The scalar or pseudo-scalar field 
4> is assumed to be a spatially homogeneous and isotropic background for the EM field. In 
Coulomb gauge, d l Ai = = Ao, the Euler-Lagrange equations then become 

f 2 A'[ - f 2 V 2 Ai + 2ff'A' i + g'e^djAk + a 2 m 2 Ai + dA\ = , (3.4) 

where €123 = +1- Notice that in the case of a non-trivial f{(fi) the value of the conductivity 
can or can not have an explicit factor of f 2 as well, depending on whether the factor f 2 
embraces the entire Lagrangian including the photon vertex or only the kinetic term F 2 as 
in eq. (3.2). 

Before analysing the impact of the modifications to EM at late times, it is instructive 
to look at the standard mode equation for f(4>) = 1, g(<j)) = = m(0). Upon transforming 
to Fourier space, the vector potential is expressed in the helicity basis, £^ with h = ± and 
we can readily write down the mode equation as 

A'l + aA' h + k 2 A h = . (3.5) 

Assuming the friction term to be time-independent for simplicity (the reasons why this as- 
sumption does not affect the general results will become clear through the section), we can 
write the solution as 
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where k = k/a. We can study the behavior of the solution in the two separate regimes, 
k 1 and k> 1, which leads to 

A h ~ Ci e~ &r i + C 2 e~ R2 ^ for 2k < 1 , (3.7) 

A ~ e"^/ 2 (Ci e- <fc " + C 2 e ife ") for 2k > 1 . (3.8) 

We recover the known result that the vector potential is either essentially frozen in the plasma 
or a rapidly oscillating and decaying wave, respectively. The transition between these two 
regimes is given by the comoving conductivity of the plasma a. If we are interested in the 
Universe around and after recombination until today, it is easy to see that the conductivity 
is tremendously large, being still a ~ 10 20 /i/Mpc at z = 1, see lower panel of figure 1. 

3.1 Running fine-structure constant 

We focus here on the case with g((f>) = m((p) = and study the effects of the term 
f 2 (cj))F a pF al3 in the action. The Fourier space equation we want to solve is thus 

fAl + &A' h + 2ff'A' h + f 2 k 2 A h = , (3.9) 
A!l l + aA! h -fjA h -a f jA h = Q, (3.10) 

where in the second line we have defined Ah = fAh and have also made the substitution a —> 
f 2 a to exemplify the scenario where the factor f 2 appears in front of the entire Lagrangian. 

In the latter case we can imagine solving this equation for /'// nearly constant (and 
negligible /"//) so that it reduces to something very similar to the previous case. This 
is realised by choosing / = /oexp{A<r(?y — rjo)} where A is a dimensionless constant, rjo is 
conformal time today, and we chose /o = /(^>(f?o)) = 1- The solutions to eq. (3.10) now read 

A h = d (5 -|{ 1 + 2A +V / ( 1 + 2A ) 2 - 4k2 }^ + Q 2 e -|{l+2A-V(l+2A) 2 -4 K 2}^ ^ ^ 

We see immediately that on large scales k <C 1/Mpc the solutions are going to behave in 
much the same way as in the case without coupling, as long as A > —1/2. Furthermore, 
if we imagine to retain only the coupling f 2 F 2 and not the f 2 j a A a then the fine-structure 
constant will be running with the scalar field. For a similar /, the running is given by 

— = -2^ « 2|AaA77| . (3.12) 
a f 

An overview on the upper bounds on the variation of a can be found in ref. [36]. In any case, 
even a conservative limit such as Aa/a < 10 -5 poses a very strong constraint on the effects 
of the coupling under inspection, i.e. |A| < 10~ 25 (Ary/Mpc) -1 . 

3.2 Running axial coupling 

This case corresponds to f(4>) = 1 and m{<f)) = 0. The equation of motion for the vector 
potential now reads 

A'h + a A h + (k 2 + hkg') A h = . (3.13) 
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Obviously, one can study this equation in great detail by numerical means. But as in the 
previous case, for an order of magnitude estimate it is sufficient to look for analytical solutions 
for roughly constant go = g' . We define £ = go /a and find 



A h = Ci e "H 1+ ^ 1_4h<_4K2 }^ + C 2 f rk{ 1 -^/ 1 - Ah '<- AK ' 2 }^ 



(3.14) 



An axial coupling of this kind gives rise to cosmic birefringence and dichroism, see for 
instance ref. [37]. Hence, the most important effect is the coherent rotation of the linear 
polarisation vector for polarised light that is traveling to us from the last scattering surface 
or from distant objects such as quasars [38]. Current observations limit such effects to be 
of the order of about a degree for CMB photons [39] or light from distant quasars [40]. In 
terms of the coupling g, the change in the polarisation rotation angle 9 is given by 

A0 oc Ag((f)(r])) = c/(r/ sourcc ) - g(rj ) ^ g' Ar] . (3.15) 

For instance, taking the CMB limit alone, the allowed average rotation is AO < 0.02 rad. With 
Ar/ ~ 14000 Mpc from the last scattering surface to today, we realise that go ^ 10 -6 Mpc -1 . 
For an average conductivity of a ~ 10 22 Mpc _1 , we have £ < 10 -28 in eq. (3.14) and any 
effect of this coupling evanesce due to the much more important conductivity friction term. 

In principle, one can imagine to get around this problem by employing some kind of 
resonance in the coupling [41]. As a representative example, let us consider the case of 
an oscillating coupling of the form g'(r]) = go sin(o;r/), where u represents the oscillation 
frequency (in conformal time). To grasp the broad features of such a choice, we can roughly 
approximate the sine with a square wave and use the solutions for constant g' as before. The 
solutions now grow periodically for one half of the cycle and decay for the other half with 
alternating cycles for opposite helicities. If, for example, there is a slight inequality of ups 
and downs one of the helicity modes will be amplified as compared to the other, on average. 
A half-cycle lasts for r/ c ~ tt/uj, and from eq. (3.14) we see that the net amplification within 
a half-cycle is exp{— hK^arj c }. For a source that is Ar] away from the observer, a photon goes 
through about n c = Arj/r] c half-cycles, and the typical difference between ups and downs is 
1/ \ffic. All together we obtain an effective exponential growth given by 

\A h \ « C 2 e~ hnA/2,i9o/ ^^ J] ^ ) . (3.16) 

Focusing on large wavelengths, for instance 1 Mpc, then k/a ~ 10~ 22 , while for CMB photons 
we still have go < 10~ 6 Mpc -1 . To make the exponent larger than unity for the negative 
helicity mode, we would need a tiny ui: 



7T (k . \ 2 / 3 

Ar/ \a 



fk V /d 

f-rffoAr/J ^lO-^Mpc" 1 . (3.17) 
The conclusion is that there is essentially no oscillation and we are back to constant g' . 



3.3 Photon mass 

In this case, with the choice /(</>) = 1 and g(4>) = 0, we allow for a non-zero photon mass m. 
Defining rh = am, the equations of motion read 

Al + aA' h + (k 2 + m 2 ) A h = (3.18) 
X " + a X ' + 2k 2 ^ X ' +(k 2 + m 2 ) X = 0- (3.19) 
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The second equation governs the evolution of the longitudinal degree of freedom x, which 
necessarily appears when the photon acquires a mass. In the simplest case of a constant mass 
term the solutions for the three polarisations coincide 

(A h , x) = c x e -§{Wi-4« 2 -vH + C2 e -i{i+vi-4~ 2 -v}*, ; {3 20) 

where \i = rh/a = mjoE- 

A recent review for constraints on the photon mass is ref . [42] , where from we learn that 
a conservative limit on rh today (which is m itself) is rh < 10 -18 eV which corresponds to 
10 u Mpc _1 . This means that \x < 10 9 today. Once again, we conclude that the conductivity 
is too large for the photon mass to have any sizable effect on the solutions. Things would 
hardly be different for a slow time- variation of rh, for the only term where this intervenes in 
eq. (3.19) is suppressed by a factor of (k/rh) 2 . 

4 Summary and conclusions 

Large-scale magnetic fields are present not only in bound cosmic structures such as galaxies 
and clusters, but also they appear to fill a large fraction of the intergalactic space. This 
fact is difficult to explain by means of standard astrophysical processes. Also the early 
Universe mechanisms of magnetogenesis are mostly not viable or not efficient, even if one 
allows for non-minimal EM. On the other hand, the non-linear dynamics of the cosmic plasma 
in the recombination era do source large-scale magnetic fields, even though with insufficient 
amplitudes. Taking this mechanism as a promising starting point, we have carefully revisited 
the non-linear dynamics in a more general setting, allowing for non-minimal couplings of EM 
to (effective) scalar fields, in the hope of boosting the magnetic seeds. 

We have derived the fully general and non- linear evolution equation of the current, and 
specialised it to the photon-electron-proton plasma, relevant in the recombination era, taking 
into account all possible effects up to second order in perturbations. We paid particular 
attention to allowing for possible effects of non- minimal couplings of EM, refraining from 
applying the Maxwell and Einstein field equations that we want to modify. Subsequently, we 
have discussed the relevance of the different effects on the current and derived a cosmological 
Ohm's law that is valid even in non-minimal EM. We also point out that Ohm's law is 
valid, and the resistivity and conductivity are well defined, on large cosmological scales, even 
beyond the Hubble horizon, as long as super-horizon correlations were initially imprinted 
into the cosmic plasma (through inflation). 

We have investigated several different types of couplings of EM to scalar and pseudo- 
scalar fields. Candidate fields which could do the job are for example the dark energy field, 
or a cosmological axion background. We have explored the possibilities of a running fine- 
structure constant, f 2 ((p)F 2 , a running axial coupling, g((p)FF, and an effective photon mass 
m 2 (4>)A 2 . In all three cases we conclude that current observational and experimental limits 
do not allow for additional amplification of EM fields on cosmological scales, because of the 
large conductivity of the plasma, even at the largest scales. 
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A Conventions 

We work with a metric signature ( h + +)• For tensor components, Greek indices take 

values ... 3, while Latin indices run from 1 to 3. We employ natural Heaviside-Lorentz 
units such that c = h = fee = eo = Po = 1 an d the electric charge is e = y/Aira. The 
reduced Planck mass is defined as Mp\ = (8-^(7) _1 / 2 . Where the metric is not specified we 
apply the 1+3 covariant Ehlers-EUis formalism [43, 44], defining a splitting of space-time 
and all tensors w.r.t. a time- like 4- vector field u^, with u 2 = u a u a = —1, referred to as a 
generic observer frame. If not mentioned explicitly, we do not specify any further. The 
covariant derivative is split into a time derivative parallel to denoted by a dot, and a 
spatial derivative, V M , projected with h^ u = g^ u + u^u u orthogonal to u M . We mostly follow 
ref. [45] for the notation, which we summarise in appendix B. 

When using a Friedmann-Lemaitre- Robertson- Walker (FLRW) metric, the Hubble para- 
meter is H = a/a where the over-dot denotes the derivative w.r.t the cosmic time t. The con- 
formal Hubble parameter is % = a' /a, where ' denotes the derivative w.r.t the conformal time 
r) defined as dt = adr), a being the scale factor normalised to unity today so that comoving 
scales become physical scales then. The Hubble parameter today is Hq = 100 /ikm/s/Mpc. 
For numerical estimates we use the ACDM model with parameters estimated from WMAP-7, 
ACT-2008, Hq an d BAO data given in table 5 of ref. [46]. The total matter and baryon dens- 
ity fractions are O m = p m /pc and f2& = Pb/pc, with the critical energy density p c = 3-Hg.Mpj. 
We also use the conversion 1 Gauss ~ 4.8 x 10 56 Mpc~ 2 in Heaviside-Lorentz units. 

One can introduce the helicity basis as 

4 = -^(«?±fe5) , (a.i) 

where (e^ , , k\ form an orthonormal spatial comoving basis with |e,^| 2 = l, k = k/k. 
In the helicity basis, the vector potential takes the form -4+£+ + A-E-. Thus, one can 
study the evolution of the Fourier modes, Ah(r], k), with respect to the helicity basis for the 
polarisation states, h = ±. 

If the generated magnetic field is statistically homogeneous and isotropic, its spectrum 
is determined by two scalar functions Ps(k) and P^k). Due to the divergence-free nature 
of the magnetic fields, the two-point function for a comoving wave vector k in Fourier space 
can be written as 

(Bi( V ,k)B;(ji,q)) = ^^5{k-q){(5 l3 -kk)P S {ri,k)-i€ i]n k n P A {^k)} (A.2) 

where Ps and Pa are the symmetric and anti-symmetric parts of the power spectrum, re- 
spectively. The symmetric part of the spectrum determines the energy density: 

(BifaQBtfaq)) = (2vr) 3 J(fc - q)P s ( V , k) . (A.3) 
With respect to the helicity basis the spectra can be directly written as 

P S /A(v,k) = k 2 (\A+(Tj,k)\ 2 ±\A-(r,,k)\ 2 ) , (A.4) 
where the upper (lower) sign corresponds to Ps(Pa), respectively. 
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B Fluid and current evolution in the covariant approach 



Let us briefly summaries the notation of the 1+3 covariant Ehlers-Ellis formalism [43, 44]. 
We will mostly follow the notation of ref. [45] who's appendix gives a concise and quite 
complete collection of definitions, relations and relevant equations. We give the covariant 
energy-momentum conservation equations and derive the general evolution equation of the 
bulk velocity field for a plasma without anisotropic stress in its rest frame. 



-1 



B.l The covariant approach 

The covariant approach employs a time-like unit 4- velocity field, with it 2 
to define a 1+3 splitting of space-time. The frame represents a family of comoving 
"observers". Then, all vectors and tensors are irreducibly decomposed w.r.t. u^, the projec- 
tion tensor = g^ v + u^u v and the projected alternating tensor e^ ua = rj^a^u^ where 
= — V~9^[^u^a^]- The covariant derivative parallel and orthogonal to are denoted 
by a dot and V M respectively. For a generic tensor Y they read 



v x y»~:.. u = h%h% ■ ■ ■ Kv n Y 
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(B.l) 
(B.2) 



The geometry as measured by the family of observers is described by the irreducible 
components of V v u^ 



= -ii^Uy + + a^v + . 



(B.3) 



These are the kinematic or geometric quantities of u^: the expansion O = V a u a , the acceler- 
ation ii^ = ^VcUp, the shear a^ u = V/^WjA, and the vorticity uj^ u = Vr^^i. The vorticity 

-A t "py aU ^_ Here the projected symmetric trace- free parts are 



vector is uj^ 
defined as 



= h£V a , Y^ u = h%Y au , 



h% h v } - \h? v h 

a p q v 



Y' 



a/3 



(B.4) 



and indices in round (square) brackets are being (anti-)symmetrised. Einstein's equations and 
the Bianchi identities then reduce to evolution equations and constraints for the geometric 
quantities, see references given above. 

Next, we decompose the energy-momentum tensor of the species s, Tg U , w.r.t. into its 
dynamic quantities or fluid variables: energy density p s = u a upT^ , isotropic pressure (inch 



possible bulk viscosity) p s 



U rpa/3 



energy flux (inch possible heat flux) g£ 



-u n T> 



(fj,)a 



and anisotropic stress (inch possible shear viscosity) ns u 
have 



In a generic frame we then 



Tr = p s u»u v + p s h^ + 2q^u u) + vr; 



IIV 



(B.5) 



even if the species is well described by a perfect fluid. The reason why q$ and lis " appear is 
that ii M is not generally the rest frame of the species s. 
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B.2 Observer frames 



Several choices of frame are possible and discussed in the literature. The physics do not 
depend on the choice of frame, but one frame may be more practical than the other. In 
the context of plasmas, two frames are interesting to note: the particle (or Eckart) frame in 
which the particle flux of the plasma vanishes, and the energy (or Landau) frame in which 
the energy flux vanishes. In the absence of heat flux, i.e. in thermal equilibrium, the two 
frames are equivalent. In this case we use the term rest frame for the frame, Us , defined by 
requiring = qt 11 = —u sa Ts^ a . For more details on these frames see for instance [27]. 

Finally, the algebraic transformations of the dynamic quantities from one frame to 
another are known exactly, see e.g. the appendix of [45]. Here we only give the transformations 
from the rest frame to a generic frame, for a species without heat flux nor anisotropic stress 
in its energy frame. One decomposes the rest frame w.r.t. the generic frame according to 

u{f = 7s (u" + t£) , V ? = , ls = (l-v 2 s y 1/2 , (B.6) 

where 7 S is the local Lorentz factor of the transformation, and Vg is the relative velocity. 
Then the dynamic quantities measured in the generic frame are 

Ps = P*+l 2 s v 2 s { P * s +p* s ), (B.7) 

Ps=P*s + H V siP*s+P*s)> ( B - 8 ) 

<tf = ls(P*s+P*sK, (B.9) 

*r = t?(p:+p:)*M, (b.io) 

and remember we assumed, q* s ^ = = 7iy , in the rest frame. The other way around it can 
be seen that 

* 2 Ps + Ps m 1 1 \ 
Ps = Ps~v s x , (B.ll) 

* 1 2 Ps + Ps /D 1n \ 

P s = Ps ~ ^v s x , (B.12) 
3 1+3^ 



and the equalities 



are useful as well. 



2 / * i *\ i * Ps ~\~ Ps CD -to\ 

IsiPs+Ps) = Ps+Ps = ,,10 ' ( R13 ) 

1 + ±vj 



B.3 Energy-momentum conservation and the velocity evolution equation 

The energy- momentum conservation for a certain species s is written as V a Ts a = Cg ■ As a 
consequence of the Bianchi identities the total energy-momentum tensor is covariantly con- 
served and therefore Cs = 0. Projecting the conservation equation parallel and orthogonal 
to we have 

p s + @(ps+ Ps) + (2ii a + V tt ) Qs + aapnf = C s , (B.14) 
# + + (> Q + a^ a )<£ + (p s +p s )u» + V^ s + (u a + V Q )vrr = C<"> , (B.15) 
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where C s = —u a Cf and = haCf are the energy and momentum transfer rates from all 
other species to species s. The momentum conservation (B.15) is an equation for the energy 
flux, q% , measured in the generic frame u^. To derive the evolution equation for the velocity, 
Vs, of the s species w.r.t. it^, we use the transformations (B.7) - (B.10) and the energy 



conservation (B.14). We notice that q% 
1 



(p s +p*)vs- Then we arrive at 1 



» 



3 



(l-v 8 )Q+[u a + a a/3 









)vf-vfV a 


- 





i+M 

Ps+Ps 



V»p* a 



(B.16) 



Notice that p* is the pressure as measured in the rest frame of the s species and is given in 
eq. (B.12) in terms of the dynamic quantities in the generic frame . See also section B.5 
on pressure perturbations. 

B.4 Current evolution equation 

Again, we follow ref. [27] . The particle number flux of each species s in the generic frame 
decomposes as 



N? = n s K + «£ 



(B.17) 



with n s being the number density as seen in it'*. Thus, particles with charge e s give rise to a 
partial charge current of js = e s Ng that we decompose as follows jjf = Q s u^ + Js with the 
partial charge density and spatial currents being 



Qs 



e s n s 



egrisv 1 : 



(B.18) 



The total charge density and spatial current are Q = Q s and J M = ^ s Jf. 

To get the evolution equation of the total current, we take the time derivative of the 
spatial current 



h s vt + n s v 



</■'■> 



(B.19) 



The evolution equation for the velocity is given above, but we also need an equation for the 
number density. Let us write the particle number conservation is as V a N" = n s A s , where 
A s is the rate of s particles added or removed due to inelastic processes. By projecting the 
covariant derivative we find 



1l s 



n.s- 



vfV n n x 



(B.20) 



1 Nearly the same expression was found in ref. [27], see their eq. (33). The differences are of order vl and 
are caused by incorrect frame transformations, their eqs. (11) and (12). Even though the terms add up to the 
correct energy-momentum tensor in their eq. (9), the anisotropic stress in their eq. (12) is not trace-free as 
it should. The trace needs to be subtracted and added to the pressure, see our (B.10) and [45]. The authors 
of ref. [27] notice this in the text and carry on, but later do not realise that this leads to a certain confusion 
about the pressure when using the energy-momentum conservation equations, see their eqs. (25) and (26). 
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Finally, we insert eqs. (B.20) and (B.16) into (B.19) to find 



+ 



4 
3 


— 


?{ 


"1 

-V 

3 


l + \v 


2 

5 


Ps+Ps L 


1 + lv 


2 

5 


Ps+Ps L 



+ (T a pvfvt 



J?-V a (v?J?) 



e s n s V»p* s + p* s J» 



(B.21) 



The first line describes geometrical effects on the total current: the first and second terms are 
due to isotropic, shear and vortical expansion, while the third term is due to the local charge 
separation coupling to the acceleration of the frame u^. The second line describes non-linear 
dynamical effects from the multi-component plasma. The third line describes the effects from 
pressure fluctuations that are known as the Biermann battery. The last line is due to the 
interactions of the plasma. The first and second term are the momentum and energy transfer 
between the species, respectively. They also contain the electromagnetic forces that give rise 
to the Hall effect and ultimately to the generalised Ohm's law. The third term comes from 
particle creation and annihilation adding or removing charges. 

B.5 Pressure perturbations 

To model the pressure in the rest frame we define two relevant quantities: the temporal and 



the spatial sound speeds, ct and c s , respectively: 

< V a p s = c? s u a V aPs , V£p* = c 2 s 



(B.22) 



Here it°V a is the time derivative in the rest frame Us and Vs is the spatial derivative 
projected orthogonal to its . Both sound speeds describe how the pressure reacts to changes 
in the density. In addition, the temporal sound speed encodes possible changes in the pressure 
due to the time evolution of the entropy, while the spatial sound speed also contains the effect 



of spatial entropy fluctuations. In case of adiabatic time evolution, q 



with c a being 



the standard adiabatic sound speed. In the absence of spatial entropy fluctuations, c s = c a , 
as well. 

We are interested in the combination V^p* + p* s Vg that appears in the velocity and 
current evolution equations above. We transform the derivatives, introduce the sound speeds 
and transform everything back into a generic frame vf 1 . We find 



^s,s 



1 + c 



) V" Ps + Ps 



Ps+Ps 



v"W) + (t£ji 



clAv»p s -(Ps+Ps) e< + (V a <)< + V^ s 2 ) +C c < + 0(3) (B.23) 



where, in the second step, we use the energy conservation equation (B.14) and took along 
only terms up to second order in perturbations. 
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